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ABSTRACT 

We present an analysis of gamma-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi 
Gamma-ray Space Telescope in the region around SNR S147 (Gl 80. 0-1. 7). A spatially extended gamma-ray 
source detected in an energy range of 0.2-10 GeV is found to coincide with SNR SI 47. We confirm its spatial 
extension at >5cr confidence level. The gamma-ray flux is (3.8 ± 0.6) x 10~^ photons cm~^ s~^ corresponding 
to a luminosity of 1.3 x 10^"^ (J/ 1.3 kpc)^ erg s~^ in this energy range. The gamma-ray emission exhibits a 
possible spatial correlation with prominent Ha filaments of SI 47. There is no indication that the gamma-ray 
emission comes from the associated pulsar PSR J0538+2817. The gamma-ray spectrum integrated over the 
remnant is likely dominated by the decay of neutral tt mesons produced through the proton-proton collisions 
in the filaments. Reacceleration of pre-existing CRs and subsequent adiabatic compression in the filaments is 
sufficient to provide the required energy density of high-energy protons. 

Subject headings: acceleration of particles — ISM: individuals 147) — radiation mechanisms: non-thermal 



1. INTRODUCTION 

GeV gamma-ray sources associated with middle-aged su- 
pernova remnants (SNRs) interacting with molecular clouds 
have recently been discovered with the Large Area Tele- 
scope (LAT) onboard the Ferm i G amma-ray Space Tele - 
scope (lAbdo et al.l 120091 l20ToM : ICastro & Slan"a l20Toh . 
The GeV emission of these SNRs is plausibly dominated by 
the decay of tt^ mesons created by proton-proton collisions, 
although a dominant electr on bremsstra hlung component is 
an alternative interpretation (lAbdo et al.ll2009i) . The SNRs in- 
teracting with molecular clouds are found to be more lumi- 
nous gamma-ray sources (~ io^^~^^ erg s ~^ in the L AT band) 
than other types of SNRs (see Thompson e t al. 1120 11) . The ob- 
served luminosity requires high-density gas so that collisions 
between relativistic particles and target gas become efficient 
enough, which is easily explained by interactions with molec- 
ular clouds. Two different types of scenarios are put forward 
to explain the bright gamma-ray emission. One scenario ex- 
plains the gamma rays as the 7r^-decay emission due to in- 
teractions between nearby molecular c louds and relativistic 
protons escaping from the SNR system (lAharonian & AtoyanI 
|T996). Another scenario considers the tt^ -decay emission 
from shock-compressed clouds in a middle-aged SNR, where 
the accelerated particles frozen in the shocked clouds effi- 
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ciently collide with target gas (lUchiyama et al.ll20TQl) . 

Recently, GeV gamma-ray emission from the Cygnus 
Loop, a middle-aged remnant without a clear signature of in- 
teractions with molecular clouds , has been detected with the 
Fermi LAT (Katagiri et al]|20li1) . The gamma-ray luminos- 
ity is quite low (^ lO"^-^ erg s~^) compared to the SNRs with 
molecular cloud interactions. The gamma-ray emission can 
be interpreted to come from diffuse gas behind the blast wave, 
but dense optical filaments seen in an Ha image would also 
offer a plausible site of the gamma-ray production. Detections 
of other low-luminosity SNRs with the Fermi LAT, thanks to 
an ever increasing exposure time, will help understand the 
gamma-ray production sites and constrain cosmic-ray (CR) 
acceleration in such SNRs. 

SNR S147 (Gl 80.0-1. 7), located toward the Galactic an- 
ticenter, is one of the most evolved SNRs in our Galaxy. 
SI 47 has a nearly circular shape with an angular diame- 
ter of 200'. It contains a complex network of long fil- 
aments in the optical band (iLozinskayal |1976'; Drew et al. 
2005), which are also visible in the radio band (Kundu et al. 
119801: ISofue et al.1 119801: IXiao et al.1 l2008h . These observa- 
tions also show that the radio and Ha emissions correlate 
very well. The velocity of th e optical filaments is estimated 
to be 80-120 km s ~^ (■Lozinskavalll976l: iKirshner & Arnold! 
1979; Phillips et al.l[T98lir No indication of interactions with 
molecular clouds has been reported. 

The radio flux density is 70 Jy at 1 GHz, and the radio 
spectrum i s known to have a spectral break around 1.5 GHz. 
IXiao et all (|2008i) obtained the integrated spectrum of SI 47 
in the 0.1-5 GHz range, and determined the spectral indices 
a ^ 0.30 ±0.15 and a ~ 1.20 ±0.3 be low and above the 
spectral break, respectively. Xia o et all (12008) also found 
that the filamentary and diffuse emission in SI 47 have dif- 
ferent spectral indices of a ~ 0.35 and ~ 1.35 respectively in 
the frequencies of 2.6-4.8 GHz (above the break). The ra- 
dio emission is considered to be synchrotron radiation from 
high-energy electrons, while the origin of the spectral break 
is uncertain. No X-ray emission has been reported to date 
from this region (Sauyageot et al. 199 0) nor any TeV emis- 
sion. The EGRET detected a GeV gamma-ray source (3EG 
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J0542+2610) in the vicinity of SI 47, but its association with 
the SNR was unclear due to it s large positional unc ertainty 
(0.°70 at 95% confidence level) ("Hartman etJ]!!!!!). 

PSR J0538+2817 plausibly associate d with S147 is lo- 
cated near the center of the SNR (iNg et al. ' 2007). The 
pulsar has a spin period of 143 ms and a spin-down lu- 
minosity of 5 X lO-^"^ erg s~^ It was first discovered in 
the radio (Anderson et al." "1996) and was later also found 
in the X-rays (McGowan et al. 2003). The X-ray observa- 
tion also revealed extended emission indicative of a pulsar 
wind nebula (PWN) ( Romani & Ng 2003). Estimated pul- 
sar ages differ significant ly bet ween a kinematic age of ~ 
3 X 10"^ yr (iKramer et aH 12 003") and a characteristic age of 
6 X 10^ yr (lAnderson et al.i ri996). The kinematic age of the 
pulsar is broadly consistent with estimate s of the SNR age 
which are in the r ange of (2-10) x 10"^ yr (iSofue et al.lll98Ql: 
iKundu et al.l[T98Qh . 

The distance to the SNR is likely to be J = 1.3 kpc given 
a plausible association with pulsar PSR 105 38+28 17. The 
distance to the pul sar is estimated to be 1.30;^q kpc from 
the parallax ( C hatter] ee et al. I I2009h and 1.2 ± 6.2 kpc with 
a dispersion measure (ICordes & Lazioll2002h . On the con- 
trary, J < 0.88 kpc is suggested by absorption lines of the 
Ble s tar HD3 6665 whic h orig inates in gas assoc iated with 
S147 ("Phillips et al.l[T981i: iSalli^n & Welshl[2004l) . The di- 
ameter of the SNR is Z) 76 (J/1.3 kpc) pc. 

Here we report the LAT observations of SNR S147. A GeV 
gamma-ray source spatially coincident with SI 47 is desig- 
nated as IFGL J0538.6+2717 in the Fermi LAT First Source 
Catalog (IFGL catalog) pu blished by the Fermi LAT Collab- 
oration (lAbdo et al.ll2010dh . In this paper, we present detailed 
analysis of this LAT source with much longer accumulation 
time of about 31 months. This paper is organized as follows. 
In section 2, the observations with the Fermi LAT and the data 
reduction are summarized. Analysis of the LAT source in the 
direction of SNR S 147 is reported in section 3, establishing an 
association between the gamma-ray source and SNR S 147. In 
section 4, we present modeling of the gamma-ray emission 
coming from SNR SI 47 and discuss its implications to the 
CR acceleration in middle-aged SNRs. 

2. OBSERVATION AND DATA REDUCTION 

The Fermi Gamma-ray Space Telescope was launched on 
2008 June 1 1 . The LAT onboard Fermi is a pair conversion 
telescope, equipped with solid state silicon trackers and ce- 
sium iodide calorimeters, sensitive to photons in a broad en- 
ergy band from 0.02 to >300 GeV. The LAT has a large ef- 
fective area (~8000 cm^ above 1 GeV for on-axis events), 
instantaneously viewing '^2.4 sr of the sky with a good angu- 
lar resolution (68% containment radius better than ~1° above 
1 GeV). Details of the LAT instru ment and data reduction are 
described in Atwo od et al.l (120091) . 

The LAT data used here were collected for about 3 1 months 
from 2008 August 4 to 201 1 March 1 . The Diffuse event class 
was chosen and photons beyond the earth zenith angle of 105° 
were excluded to reduce the background from the Earth limb. 
The P6_V1 10 instrument response functions were used for 
the analyses in this paper. 

3. ANALYSIS AND RESULTS 

^ http ://w w w. slac . Stanford . edu/ exp/ glast/ groups/canda/ archive/ 

pass6vl l/lat_Perf ormance.htm 



We utilized gtlike in the Science Tools analysis packag43 
for spectral fits. With gtlike, a binned maximum likelihood 
fit is performed on the spatial and spectral distributions of ob- 
served gamma rays to optimize spectral parameters of the in- 
put model taking into account the energy dependence of the 
point- spread function (PSF). 

Analysis using gtlike is performed in a 17° x 17° region 
around SNR SI 47, which is referred to as a region of in- 
terest (ROI). The fitting model includes other point sources 
whose positions are given in the IFGL catalog Since the 
nearby SNR IC443 is spatially resolved by the LAT, we model 
its spatial distribution as a disk centered at (a, ^)=(94.°31, 
22.°58) with a diameter of 48' according to lAbdo et al.] 
(1201 Obi) . The Galactic diffuse emission is modeled by 
"gll_iem_v02_P6_Vll_DIFFUSE.fit" and an isotropic com- 
ponent (instrumental and extragalactic diffuse backgrounds) 
by "isotropic_iem_v02_P6_Vll_DIFFUSE.txt". Both back- 
ground models are the standard diffuse emission models re- 
leased by the LAT tearrE|. In each gtlike run, all point sources 
within the ROI and diffuse components in the model are fit- 
ted with the normalization left free. The spectral shapes are 
either fixed or set free depending upon the specific analysis 
(see below). Note that we also include IFGL sources out- 
side the ROI but within 14° from SI 47 with their parame- 
ters being fixed at those of the IFGL catalog. The count map 
of the ROI above 1 GeV (Figure □ left) shows that the Crab 
Nebula and SNR IC443 are the dominant gamma-ray sources 
in this region. We model the Crab spectrum using three 
components (the Crab pulsar, the Crab PWN synchrotron, 
and the Crab PWN inverse Compton components), and fix 
their s pectral shapes in the fit, following the previous LAT 
study (lAbdo et al.ll2010e|) . For I C443, we niodel the emission 
as a broken power law following lAbdo et aH (1201 Ob|) . and also 
fix its spectral shape. 

3.1. Spatial Distribution 

In Figure [T] we show a smoothed count map of the 
ROI above 1 GeV, a corresponding background model map, 
and the background- subtracted count map. The background 
model map includes contributions from the Galactic diffuse 
emission, isotropic diffuse background, and nearby discrete 
sources. The model parameters of diffuse components and 
nearby sources are set at the best-fit values obtained by gtlike 
using the data above 1 GeV. An extended gamma-ray source 
associated with S147 is visible in the background- subtracted 
map. In the model, the spatial distribution of the S147 source 
is assumed to be the Ha image (lFinkbeineri2003|) (see below). 
The spectral shape is fitted as a power-law function with the 
index set free. The point source in the inset is simulated with 
the same spectral shape as SI 47 obtained by maximum likeli- 
hood fits. 

We generated a Test Statistic (TS) map using the LAT data 
above 1 GeV in the SI 47 region (Figure |2] left). TS is defined 
as -2 A In(likelihood) obtained by gtlike between models of 
the null hypothesis and an alternative. In this paper, we refer 
to TS as a comparison between models without a target source 
(null hypothesis) and with the source (alternative hypothesis) 
unless otherwise mentioned. In this map, the TS value at 
each grid position is calculated by using a model with a point 

^ available at the Fermi Science Support Center. 
http://fermi.gsfc.nasa.gov/ssc | 

^ The flux of SI 47 would change < 10% below 1 GeV and < 4% above 
1 GeV when we use 2FGL catalog ( Nolan et al. 2012). 



Gamma-ray Emission from SNR SI 47 



3 




Fig. 1 . — (Left) Fermi-LAI count map above 1 GeV around SNR SI 47 in units of counts per pixel. The pixel size is OP 1. Smoothing with a Gaussian kernel of 
a = 0?25 is applied. SNR S147 is represented by a white circle. The background sources contained in the ROI are shown as green crosses. (Middle) Background 
model map. The green boxes each with the dimensions oflPOx 1?5 represent the regions used for the evaluation of the accuracy of the Galactic diffuse model. 
(Right) Background- subtracted count map. A simulated point source is shown in the inset. 




Fig. 2. — Left: TS map obtained with maximum likelihood analysis of the Fermi-LAI data in the vicinity of SNR S147 above 1 GeV. Overlaid are linear 
contours of the background- subtracted count map above 1 GeV taken from Figure [T] A white circle represents the outer boundary of SNR S147. A cyan cross 
and circle represent a position and positional error (95% confidence level) of 3EG J0542-I-2610, respectively (Har tman et al.iil99^) . A magenta cross indicates 
the position of PSR J0538-1-2817. Right: Ho; flux intensity map of SNR S147 in units of rayleighs (10^/47r) photons cm~^ s~^ sr"^ cFinkbeiner 2003) , with the 
contours of the background- subtracted count map overlaid. 
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Fig. 3. — Five divided regions for comparison of gamma-ray and Ha emission. The grids overlay (Left) the LAT count map i n the vicinity of 
SNR SI 47 above 1 GeV with the background subtracted, and (Right) an Ha image convolved with the LAT PSF (see Section [3Jl . The count 
map and image is in units of counts per pixel with the size of 0.°05. Here the LAT count map is smoothed by a Gaussian kernel of a = 0.°25, to 
show the clear spatial distribution of the S147 source with the low number of photon counts. Note that the the count map is not smoothed for 
the comparison of each region. 



source placed at the position, which has a power-law energy 
distribution with its index being free. The excess gamma rays 
above backgrounds are distributed inside the SNR boundary 
and the spatial extent is consistent with the remnant size. To 
evaluate the spatial extent, we fit the data with a model where 
the spatial distribution of the SI 47 source is set to a disk im- 
age with a uniform surface brightness, the size of which cor- 
responds to that of the SNR in the Ha image (the white circle 
in Figure O. Indeed, TS obtained with this model is larger by 
~ 50 than the largest TS of a point source in the S147 region 
above 1 GeV, which means the SI 47 source is extended with 
>5a confidence. The spectral distribution of the S147 source 
is assumed to be a power-law function with its index free. 

In Figure [2] (right), the background- subtracted gamma-ray 
map is compared with an Ha emission map of SI 47. To per- 
form a more detailed comparison between the gamma-ray and 
Ha maps, we divide the S147 region into five cells as shown 
in FigureO and compare gamma-ray (> 1 GeV) and Ha fluxes 
in each cell. We note that the size of the cell (~ 1°) is com- 
parable to the LAT angular resolution (better than ~ 1 ° above 
1 GeV). To take into account the LAT PSF for the morpho- 
logical comparison, the Ha image is convolved with the LAT 
PSF as shown in Figured The gamma-ray map used for this 
comparison is a count map with all model components except 
for the S147 source subtracted as described above. Since the 
PSF depends on energy, we calculate the PSF assuming the 
spectral index of the SI 47 source to be 2.5, the best-fit value 
obtained with the Ha image used as a template above 1 GeV. 
Note that the exposure of the LAT observations is uniform 
within 0.5% in this region. 

Figured] shows the resulting correlation diagram. The plot- 
ted gamma-ray counts are obtained by summing up the sub- 
tracted count map of the LAT in each region, while the Ha 
counts are calculated by summing up the Ha image convolved 
with the LAT PSF. The figure shows a possible correlation be- 
tween gamma-ray and Ha fluxes. Though it is challenging to 
establish the presence of the correlation given large statistical 



uncertainties, the diagram suggests that an extension of the 
Fermi mission (~ 10 yr) will provide an opportunity to con- 
firm it. 




normalized Ha counts 



Fig. 4.— Relation of the LAT (1-200 GeV) and Ua counts. Note 
that both are normalized to make their total counts 1.0. The blue line 
represents a linear function with its slope of 1.0. The error bars repre- 
sent statistical uncertainty associated with the LAT-observed photon 
counts (Icr). 

We also test different spatial templates, a disk template, 
sphere template, and shell template, as shown in Figure|5] The 
disk template is a disk with uniform surface brightness. The 
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sphere template is a two-dimensional projection of a sphere 
with a uniform emissivity per unit volume. The shell is a 
two-dimensional projection of a spherical shell with a ratio 
of inner diameter to outer diameter set to 0.8 based on the 
Ha map. Center positions and diameters of the disk, sphere, 
and shell templates are fitted to the data. Using the different 
templates, we perform maximum likelihood fits and compare 
the best-fit parameters in the energy range of 1-200 GeV. The 
spectral shape of the SI 47 source in the wide energy range is 
assumed to be a power-law function. As shown in Table [T] 
the SI 47 source is significantly detected in each case, and the 
obtained fluxes and spectral shapes are almost the same. The 
Ha image has the largest TS among all templates, despite the 
fact that the other templates have three more free parameters 
(position and diameter) than those of the Ha image. 

Given the results of the correlation diagram and the fits of 
the different spatial templates, the gamma-ray emission from 
the S147 source has a possible spatial correlation with the Ha 
filaments. Hereafter we adopt the Ha image as the spatial 
template of the SI 47 source. 

3.2. Spectrum 

The spectral energy distribution (SED) of the source asso- 
ciated with S147 is obtained by dividing the 0.2-200 GeV 
energy band into six energy bins. Since each energy range for 
the fit is small, we model the Crab as a point source with a sin- 
gle spectral component instead of three. To model the bright 
emissions from Crab and IC443 in the small energy range, 
these spectral shapes are adopted to be power-law functions 
with their indices free. The indices of the other sources are 
fixed at the values in the IFGL catalog. S 147 is fitted with a 
simple power-law function in each energy bin with its photon 
index fixed at 2.1, which is obtained by a broadband fitting in 
0.2-200 GeV (see below). We note that the flux obtained for 
S 147 in each energy bin is insensitive to the choice of index, 
if it is confined to a reasonable range (say, 1-3). In each fit, all 
sources within the ROI and diffuse components in the model 
are fitted with their normalization being free. 

The systematic errors in the spectral analysis are mainly 
due to uncertainties associated with the underlying Galac- 
tic diffuse emission and uncertainties of the effective area of 
the LAT. The uncertainties of the Galactic diffuse emission 
are primarily due to the imperfection of the Galactic diffuse 
model and/or the contributions from discrete sources not re- 
solved from background. We evaluate the uncertainties of 
the Galactic diffuse emission by measuring the dispersion of 
the fractional residuals in 10 regions around SI 47, where the 
Galactic diffuse component dominates (Figure [T]). The frac- 
tional residuals, namely (observed-model)/model, are calcu- 
lated in three energy bands (0.20-0.45 GeV, 0.45-1.0 GeV, 
and 1.0-10 GeV) for each region. At each energy range, 
the uncertainties of the Galactic diffuse model are adopted 
separately as the second largest values among 10 residuals 
(90% containment). From the results, the uncertainties are 
evaluated as 5.3%, 6.4%, and 8.9% at energy ranges of 0.20- 
0.45 GeV, 0.45-1.0 GeV, and 1.0-10 GeV, respectively. Sys- 
tematic uncertainties of the effective area are 10% at 100 
MeV, decreasing to 5% at 560 M eV, and increasing to 20% 
at 10 GeV and above (Rando et al. 2009). Note that system- 
atic errors associated with choice of spatial models of SI 47 
are negligible, since the best-fit parameters obtained by gtlike 
with the different spatial templates for SI 47 are almost the 
same as described in Section [3TT] 

Figure [6] shows the resulting SED for SI 47. Total system- 



atic errors are set by adding in quadrature the uncertainties 
due to the Galactic diffuse model and the effective area. Note 
that upper limits are calculated for the spectrum of SI 47 be- 
low 1 GeV given the conservative systematic errors. Inspec- 
tion of the figure indicates a spectrum steepening above a 
few GeV. To evaluate the significance of the steepening, we 
performed likelihood-ratio tests between a power-law func- 
tion (the null hypothesis), and either an exponentially cut-off 
power-law or a smoothly broken power-law function (the al- 
ternative hypotheses) for 0.2-200 GeV data. The exponen- 
tially cut-off power-law function is described as 



dN 
dE 



■■No 



exp 



E 



(1) 



where Eq is 1 GeV. The photon index F, a cutoff energy 
Ecuu and a normalization factor No are free parameters. The 
smoothly broken power-law function is described as 



dE \Eo 




r2-ri 



(2) 



where Eq is 1 GeV. The photon indices Fi below the break, F2 
above the break, a break energy £'br, and a normalization fac- 
tor No are free parameters. The resulting TS are TS cutoff = 
-21n(LpL/Lcutoff) = 2.4 and TSbpl = -21n(LpL/LBPL) = 6.8, 
corresponding to a significance of 1 .5cr and 2.1a respectively. 
There is an indication of the spectral steepening above a few 
GeV, but the possibility of a simple power law spectrum can- 
not be conclusively rejected. The parameters obtained with 
the broken power law model are photon indices Fi = 1 .4 ib 0.5, 
F2 = 2.5 ± 0. 15, and ^break = 1 .0 ± 0.8 GeV, with an integrated 
flux in 0.2-200 GeV of (3.8 ± 0.6) x 10"^ photon cm'^ s'^. 
The photon index obtained with the simple power law is 
2.14 ±0.05. The gamma-ray luminosity in 0.2-200 GeV is 
calculated as 1 .3 x 10^"^ (J/ 1 .3 kpc)^ erg s~^ using the result of 
the broken power law. The best-fit functions are represented 
in Figure [6l 

3.3. Upper Limit on PSR J0538+281 7 

As shown in Figure [2l no gamma-ray counterpart is found 
at the position of PSR J0538+2817 located near the center 
of the SNR. An upper limit to the GeV flux from the pul- 
sar is determined by performing a maximum likelihood fit in 
0.2-200 GeV. The model is the same model as in Section [T2l 
where SI 47 is modeled as a smoothly broken power-law func- 
tion with the normalization, the energy break, and the photon 
indices below and above the break set free, except for adding 
a point source at a position of PS R J0538+ 2817, specifically, 
(a, (5)=(84.°604, 28.°286) (Anders^eLanri996). The spec- 
tral shape of the pulsar is assumed to be an exponentially 
cut-off power-law function (see Equation [TJ with a photon 
index of 1.3 and a cutoff energy of 1.9 GeV (iRomani et al.l 
|2011|) . Note that these assumed parameter s are within typica l 
values of the observed gamma-ray pulsars (lAbdo et al.ll2011h . 
The calculated upper limit on the flux (0.2-200 GeV) is 1 .9 x 
10~^ photons cm~^ s~^ at 90% confidence level, corresponding 
to a luminosity limit of 4.7 x 10^^ (J/1.3 kpc)^ erg s~^ 

4. DISCUSSION 

4.1. Modeling of Multiwavelength Spectra 

We found extended GeV gamma-ray emission spatially co- 
incident with SNR SI 47 using 31 months of data acquired 
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TABLE 1 

Best- FIT Values of Different Templates for the SI 47 Source 



Template 


Center position 

(a, S) 


Diameter 


Flux 

[10"^ photons cm~^ s~^] 


Photon index 


TS 


Disk 

Sphere 

Shell 

Hq; image 


(84.°97, 27.°97) 
(85.°00, 27.°93) 
(85?29, 27?80) 


3P5 
3.°9 
3?1 


9.2±1.1 
9.2±1.1 
8.6±1.0 
8.5 ±1.0 


2.45 ±0.14 
2.45 ±0.14 
2.51 ±0.15 
2.48 ±0.14 


86.6 
84.7 
88.0 
94.8 



Note. — Flux of the S147 source is calculated in the energy range of 1-200 GeV. 




87° 86° 85° 84° 83° 87° 86° 85° 84° 83° 




87° 86° 85° 84° 83° 87° 86° 85° 84° 83° 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 5. — The S147-source templates for the LAT spectral analysis. The templates are (a) a disk template, (b) a sphere template, (c) a shell 
template, and (d) an Ha image, overlaid with the contours of the LAT background- subtracted count map above 1 GeV (the same contours as 
that overlaid in Figure |2j. 



with the Fermi LAT. The size of the gamma-ray emitting re- 
gion is consistent with the shell of Ha and radio emission 
(~2000. Moreover, the gamma-ray emission exhibits a pos- 
sible correlation with Ha filaments and therefore also with 
synchrotron radio filame nts given a tight correlation between 
the Ha and radio maps (IXiao et al.l l2QQ8). We did not find a 
gamma-ray source at the position of PSR J0538+2817, which 
is considered to be the stellar remnant of the supernova explo- 
sion. 

Let us consider the gamma-ray emissions from radio- 
emitting region s, which can be de composed into filaments and 
diffuse regions (IXiao et al.l2QQ8|) . The diffuse and filamentary 
components explain the radio data below and above a break 
around L5 GHz (~ 10""^ eV), respectively. The diffuse radio 



emission is supposed to come from the shocked gas behind the 
blast wave which is propagating in intercloud medium, while 
the filaments are formed through radiative shocks driven in 
atomic clouds. We calculate the gamma-ray spectrum of each 
zone expected from the radio spectra under some reasonable 
assumptions. These components are separated only in spec- 
tral space at gamma-ray energies, since it is difficult to spa- 
tially resolve these components with the PSF of the LAT. The 
distance to the remnant is assumed to be J = L3 kpc, which 
is the most likely distance to the pulsar PSR J0538+2817 
that is supposed to be associated with S1 47 (Ng et al. 2007; 
iChatteriee eTani2QQ9l:lKramer et al.ll2003h . Following this as- 
sumption, the radius of the remnant is 38 pc. We also adopt 
a remnant age of = 3 x lO'^ yr from the kinematic age of 
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Fig. 6. — SED of S 147 measured by the Fermi LAT. The arrows rep- 
resent upper Umit on the fluxes at 90% confidence level. The total 
systematic errors are indicated by black error bars while statistical er- 
rors (1 cr) are indicated by red error bars. The black points and arrows 
represent upper limits taking the systematic errors into consideration. 
The blue dashed, dotted, solid lines represent the best-fit power-law, 
exponentially cut-off power-law, smoothly broken power-law func- 
tions from a binned likelihood fit in 0.2-200 GeV. 

PSR J0538-F2817. The distance, radius, and age have uncer- 
tainties of 20-30% (as described in Section 1). 

After supernova explosions, when the swept-up gas be- 
comes comparable to the ejecta mass and the blast wave 
slows, SNRs enter the S edov-Taylor (adiabatic) phase (e.g., 
iTruelove & McKejll99 9). The total energy of SNRs remains 
almost constant during the Sedov-Taylor phase, because ther- 
mal and synchrotron radiation energy losses are negligible 
compared to the internal energy. The radiative phase begins 
when radiative cooling dominates the energy loss of SNRs 
and the adiabatic approximation breaks down. The transition 
time and radius to the radiative phase are 



1.3 X 10^4/ Vm yr, 



Rtr = 14 X £5^ «,CM PC, 



(3) 
(4) 



where £^51 is the explosion kinetic energy in units of 
10^^ erg and ^icm is the density of the intercloud gas in 
units of cm-^ (Ci offi et al.| [l988). If S147 is in the Sedov 
phase, the density of the intercloud gas should be ^icm < 

O.IAeK^^qI^^ cm~^ (= Her) derived from the condition of 

to < hr, where ^sokyr = ^o/(3 x 10"^ yr). This result is consis- 
tent with the density of the intercloud gas calculated by the 
relation in the Sedov-Taylor stage: 



^iCM = 3.3x 10-2^51^; 



30kyr^38pc 



cm 



(5) 



where Rsspc = ^/(38 pc). On the other hand, if we assume that 
SI 47 is in the radiative phase, the density of the intercloud 
gas should be ^icm > ^cr- This leads to a smaller remnant, 
incompatible with our preferred distance of 1.3 kpc. In this 
paper, we assume S147 is in the Sedov phase. 

The temporal evolution of the particle (protons/electrons) 
momentum distribution in diffuse and filamentary zones can 
be described by: 



|d,f 



(6) 



where be^p = -dp/dt is the momentum loss rate, and Qe§(p) is 
the particle injection rate. The super scripts "d" and "f" of the 
parameters represent the diffuse and the filamentary regions, 
respectively. To reproduce the spectral steepening in the radio 
band, we adopt phenomenological forms for injection distri- 
butions of protons and electrons: 




(7) 



where po is 1 GeV c ^ and /^^ut is a cutoff momentum. Here 
we assume the injection rate is time independent for simplic- 
ity. To obtain the radiation spectra from the remnant, N^(p^ t) 
is numerically calculated using parameters in Table |2 We 
use N^(p^to) and N^(p^to/2) to calculate the radiation spec- 
tra from the diffuse region and filaments, respectively. The 
momentum losses for electrons include synchrotron radia- 
tion, bremsstrahlung, inverse compton (IC) scattering, and 
coulomb collisions, while those for prot ons include pion pro - 
duction losses, and coulomb collisions (IS turner et al.lll997h . 
In the case of the diffuse region, the adiabatic loss is also taken 
into account using the Sedov-Taylor evolution. The gamma- 
ray emission mechanisms include the tt^ -decay gamma rays 
due to high-energy protons, and bremsstrahlung and IC scat- 
tering processes by high-energy electrons. Calculations of the 
gamm a-ray emission are done using the method described 
in Abd o et aP (120091) . The interstellar field used for cal- 
culations of IC includes an infrared blackbody component 
(^TiR = 3 X 10"^ eV, UiR = 1 eV cm"^), an optical blackbody 
component (kTo^t = 0.25 eV, Uopt = 1 eV cm~^), and the cosmic 
microwave background. The infrared and optical components 
are set to reproduce the interstellar radiation field in the GAL- 
PROP code (iPorter et al.n 2008). 

The physical parameters used for the model calculations are 
summarized in Table |2l The blast wave ve locity is determined 
by vt = OAR/to (^e.g.. iTruelove & McKeel [r999). The density 
of the intercloud gas ^icm is determined using Equation [?1 
Then, the post-shock density in the diffuse region is set by 
^Hd = 4^icM. The magnetic field in the intercloud region is 
varied within typical values of 1-5 /iG, when predicting the 
gamma-ray spectra. The compressed magnetic field in the dif- 
fuse region is determined as = \/2/ 3(nud / ^icm)^icm • 

We assume that the filaments are formed by a radia- 
tive shock wave (50 km s~^ < < 200 km s~^) driven into 
atomic clouds, which are denser than the intercloud medium. 
The post-shock gas is subject to radiative cooling, forming 
compressed gas (i.e., filaments). The gas density nuf and 
temperature Tf of the filaments are estimated fro m optical 
lines ( Kirs hner & Arnold! 119791: iFesen etal.iri985h (see Ta- 
ble [21). We ad opt Vy = 100 km s~^ according to the opti- 
cal observations (|Lozinskavalll976l:lKirshner & Amoldll 19791: 
iPhillips et al.lll98l|) . The density of the atomic clouds ^atc is 
calculated as 



^atc = ^ICM 



Vb 



= 0.83F^5i^38pcV.7"' cm-^ 



(8) 
(9) 



where Vsj = v^/(100 km s ^) and F :^ 3.2 - 4.8(v^/v/^) 
2.6(Vs/vbf (iMcKee & Cowid Il975h . The magnetic field 
strength in the filaments, Bf, can be estimated from pres- 
sure balance (.Hollenbach & McKee.1979.) . The shock ram 
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TABLE 2 

Basic Parameters of Multiwavelength 
Models 



SNR dynamics 


assumed parameters 




Distance: d 


1.3 kpc 


Age: tQ 


3 X 10^ yr 


Explosion energy: £^tot 


(1-3) xlO^i erg 


dependent parameters 




Radius: R 


38 pc 


Blast wave velocity: 


500 km s-i 


Filament gas properties 


assumed parameters 




Preshock magnetic field: 5atc 


3-10 iiG 


Gas density: ^Hf 


100-500 cm-3 


Temperature: Tf 


2 X 10^ K 


Shock velocity: 


100 km s-i 


dependent parameters 




Magnetic field: Bf 


Bf{nm,Tf,B^_ic,Vsf 


Preshock gas density: /latc 


2-6 cm~3 


Diffuse gas properties 


assumed parameters 




Preshock magnetic field: 5icm 


1-5 /iG 


dependent parameters 




Magnetic field: 


3-16 /iG 


Preshock gas density: uicm 


0.03-0.1 cm-3 


Gas density: «Hd 


0.1-0.4 cm-3 


^ see Section |4J1 



pressure Pram = ^atcMRV^^ should be balanced by magnetic or 
thermal pressure, where n^tc is the gas density of the atomic 
clouds (preshocked gas of the filaments), and /in is the mass 
per hydrogen nucleus. When magnetic pressure supports the 
filaments, the magnetic field strength is set by Pram = Bf^/^Tr. 
Using the relation of Bf = \/273(^Hf/^atc)^atc, 

Bf^Ux V,7'/'5atc-6'/'^Hf'/' /iG, (10) 

where 5atc-6 = ^atc/(l A^G). The magnetic field in the atomic 
clouds ^atc is set in a range of a few times higher than a typi- 
cal intercloud magnetic field. When the thermal pressure sup- 
ports the filaments, Pram should be equated with the thermal 
pressure XtUufkTf, where Xt = 23 assuming the filament gas is 
fully ionized and k is the Boltzmann constant. In this case, 

5f ^ 60 X V,7'5atc-67}4' /iG, (11) 

where 7f4 is 7f/(10'^ K). The lower value between Equa- 
tions [TOland[TT]determines Bf, which is found to be 5/ ~100- 
300 /iG for the adopted parameters. 

We adjust the injection spectrum of electrons Q^'^(p) (Equa- 
tion |7]) to reproduce the observed radio spectrum, which is ar- 
guably the synchrotron radiation of the relativistic electrons. 
Free parameters of Qf^(p) are af^, p^^^, and s. The index s is 
set at 1.5 to explain the observed radio synchrotron emission. 
We use the same indices for the filamentary and diffuse com- 
ponents, while the normlaization and the cutoff momentum 
are set independently. To predict the gamma-ray spectrum, 
we consider two cases where the number ratio of the relativis- 
tic electrons to protons Kep = Ge/ap are 0.01 and 1. A value 
of Kep = 0.01 is sim ilar to what is l ocally observed for CRs at 
GeV energies (e.g.. Beische r et al]l20 Q9). 

The corresponding gamma-ray spectra are predicted using 



physical parameters varied within a range described in Ta- 
ble [21 The gamma rays are dominated by the 7r^-decay emis- 
sion for Kep = 0.01. Figure |7] shows that the observed gamma- 
ray spectrum can be reproduced by using the parameters that 
are outlined above. The color-filled regions in the gamma-ray 
band represent ranges of possible fluxes with physical param- 
eters varied as described in Table |3] For case (a3), where the 
variable parameters are set at nearly central values in the ex- 
pected ranges, the total kinetic energies of protons (Wp) are 
calculated as 1.7 x 10'^'^ erg and 5.0 x 10"^^ erg for the fila- 
ments and diffuse, respectively. Although Wp of the filaments 
is much smaller than that of the diffuse, the observed gamma 
rays are dominated by the tt^ -decay emission from the fila- 
ments due to the high density. The strong magnetic field of 
210 /iG also enhances radio synchrotron emission from the 
filaments. Note that only the filament component can repro- 
duce the observed gamma-ray spectral shape in any case. 



r (al) nHf==100cm"3 




(a2) nHf=500cm-3 






(a3) nHf=250 cm-3 




1 1 



10-^ 10-5 10-3 10-1 IQi 103 10^ 10^ 10^ 10^1 



E [eV] 

Fig. 7. — SEDs of SI 47 with model curves for three cases. Cases (al), 
(a2), and (a3) represent ^Hf = 100. 500, and 250 cm~^, respectively (see Ta- 
ble [5). Kep is assumed to be 0.01. The black points represent observed data 
in the radio (Xiao et al. 2008) and LAT bands. The red and green lines rep- 
resent diffuse and filamentary components, while the blue lines represent a 
sum of the two components. The filled regions of cases (al) and (a2) in the 
gamma-ray band show ranges of the gamma-ray spectra using parameters in 
TableO Case (a3) shows sub-components of the models for a particular case. 
The synchrotron spectra from primary (dashed line) and secondary electrons 
(dot-dashed line) are drawn in the radio band, while the gamma-ray spec- 
trum can be decomposed into tt^ -decay from the relativistic protons (dashed 
line), bremsstrahlung from primary (dot-dashed line) and secondary electrons 
(vertical-dashed line), and IC scattering from primary electrons (dotted line). 



In the case of Kep = 1, the gamma-ray emission of the fil- 
aments is dominated by electron bremsstrahlung, while that 
of the diffuse region has significant contributions by both 
bremsstrahlung and IC scattering. Figure |8] shows that the 
observed gamma-ray spectrum is under-predicted for nuf = 
100 cm~^. To reproduce the observed spectrum, nuf ^ 
500 cm~^ is preferable. However, the total kinetic energy 
of the relativistic protons is only Wp = 2.5 x 10^^ erg in case 
(b3). This value is only one third of the energy content of 
Galactic CR protons in the volume of SI 47, Wqr = UqrV 
7.5 X 10"^^ erg, where Ucr ^0.7 qV cm~^ is the energy den- 
sity of CR protons and V = (47r /3)R^ is the volume of the SNR 
with a radius of 7? = 38 pc. This suggests that Kep is unlikely 
to be 1. 
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TABLES 

Parameters of Multiwavelength Models {Kep = 0.01) 







Parameters 






Energetics 




Model 




B 


Pent 




We 


Ub 




[cm-3] 




[GeVc-1] 


[10^8 erg] 


[10^6 erg] 


[eV cm-3] 


(al) filament 


100 


90-170 


60 


0.29-0.78 


0.29-0.84 


200-720 


diffuse 


0.1-0.4 


3-16 


8-15 


22-350 


42-550 


0.2-6.4 


(a2) filament 


500 


90-240 


60 


0.11-0.63 


0.080-0.56 


200-1430 


diffuse 


0.1-0.4 


3-16 


8-15 


22-350 


42-550 


0.2-6.4 


(a3) filament 


250 


210 


60 


0.17 


0.15 


1100 


diffuse 


0.2 


10 


10 


50 


88 


2.5 



Note. — The total kinetic energies (£"1011) of radiating protons (Wp) and electrons (We) are calculated for 
£'kin > 50 MeV and E'kin > 20 MeV, respectively. 



TABLE 4 

Parameters of Multiwavelength Models (Kep = 1) 





Parameters 






Energetics 




Model 


nil 


B 


Pcut 




We 


Ub 




[cm~^] 


[mG] 


[GeV c-i] 


[10^6 erg] 


[10^6 erg] 


[eV cm-3] 


(bl) filament 


100 


90-170 


40 


0.29-0.81 


0.31-0.76 


200-720 


diffuse 


0.1-0.4 


3-16 


8-15 


24-390 


46-610 


0.2-6.4 


(b2) filament 


500 


90-240 


40 


0.22-0.81 


0.18-0.76 


200-1430 


diffuse 


0.1-0.4 


3-16 


8-15 


24-390 


46-610 


0.2-6.4 


(b3) filament 


500 


90 


40 


0.81 


0.76 


200 


diffuse 


0.4 


4 


14 


250 


390 


0.4 


Note. — The total kinetic energies (£"^11) of radiating protons (Wp) and electrons (We) are calculated for 



explained by the emission from only pre-existing CR s ac cel- 
erated at the shock wave of the SNR shell. In Section l4Jl we 
modeled the radio and gamma-ray spectra without specifying 
the sources of energetic particles. In this section, we calcu- 
late the emission from the relativistic particles produced by 
acceleration (i.e., reacceleration) of the pre-existing CR s, fol- 
lowing the prescription given by lUchiyama et al.l (|2Q1Q|) . The 
emission from two regions, fila ment s and diffuse, are consid- 
ered in the same way as Section |4J1 

The spectral distributions of the reaccelerated CRs are cal- 
culated as follows. The number density of ambient CR pro- 
tons and electrons, seeds of reacceleration, are adopted as the 
spectra of the observed Galactic CR protons nGCR,p(p) and 
electrons+positrons ^GCR,e(/^)' 

nGCRAP) = '^^JpP''pf''^ (12) 

riGCRAP) = ^^JeP'-'pf (1 +/^or''', (13) 

where = p/(GcV/c), Jp = 1.9 cm~^ s~^ sr~^ GeV~\ and 
J, = 2x IQ-^ cm "^ s-i sr-i GeV-^ (IShikaze et al.1 [2QQ7I: 
IStrong etaH 1200 4). We assume that spectral distributions of 
the pre-existing relativistic particles in the atomic clouds (i.e., 
preshocked gas of the filaments) are also the same as those 
of the Galactic CRs. At the shock wave, the CRs described 
above are assumed to be accelerat ed according to the theory 
of diffusive shock acceleration (Blandford & Eichler 1987). 
The number density of the reaccelerated CRs as a function of 
momentum ^acc(/^) is described as 

riDSAip) = (a + 2) p-^ r dp'nocKip') (14) 



£kin > 50 MeV and E^^^ > 20 MeV, respectively. 



r (bl) nHf=100cm"3 

>^ , , , . 




(b2) nHf==500 cm"^ 

>^ , , , . 






r (b3) nHf=500cm"3 

'y^\ ® 


1 1 



10^ 103 
E [eV] 



109 10^1 



Fig. 8. — Same as Figure |2]but Kep is assumed as 1. Cases (bl), 
(b2), and (b3) represent nu = 100, 500, and 500 cm~^, respectively 
(see TablelU. 



4.2. Reacceleration of Galactic Cosmic Rays as Sources of 
Gamma-ray Emitting Particles 

The total energy of the relativistic electrons obtained for 
the most reasonable set of parameters (i.e., case (a3)), is 
We — 9x 10"^^ erg (see Table [3]) . This value is only three times 
higher than the total energy of the Galactic CR electrons that 
was stored in the interstellar medium SNR SI 47 has swept up. 
This indicates the possibility that the observed spectra can be 
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where ^gcr is the number density of pre- shock CRs, a = 
(rsh + 2)/(rsh- 1) and rsh is the shock compression ratio. Here 
we adopt to be 4, assuming the specific heat ratio of non- 
relativistic monoatomic gas (=5/3) and a strong shock (Mach 
number ^10). Given the finite acceleration time, the maxi- 
mum attainable momentum of the accelerated particles is also 
finite. Here we introduce a cutoff momentum (/?c_acc) for the 
spectra of the accelerated particles. Then, the number density 
of particles accelerated at the shock wave is described as 



^acc(/^) = nj)SA(p) X exp 




Since the diffuse region is behind the shock wave, the number 
density of relativistic particles in the diffuse region n^^ip) has 
the spectral distribution of Equation [15] On the other hand, 
the filaments are formed as the shocked gas is compressed 
radiatively. In this process, the particles frozen in the gas are 
heated and gain energy sls p ^ s^^^p, whe re s = (nuf/ ^atc) / ^sh , 
and th e density increases by a factor of s (|Blandford & Cowid 
Il982|) . Therefore, the number density of the accelerated and 
compressed particles in the filaments is calculated as nl p(p) = 

s^'^n,,,{s-'l^p). 

To calculate the multiwavelengt h sp ectrum, the same pro- 
cedure and parameters as Section 14.11 are applied except for 
the injection rates of the relativistic particles: 

Qip{p) = {^-f)y(n,cM/nm) x (16) 
Qipip) = fyin^ic/nm) x ^',//?)/4, (17) 

where tf^^ = to and tf^^ = to/2 are the total injection time. The 
filling factor of the preshocked atomic gas is defined as / = 
Vb/y, where Vb and V are volumes of the atomic clouds and 
the SNR, respectively. In this scenario, the relativisti c pa rticle 
distributions are almost fixed unlike those in Section 14711 Free 
parameters of the particle distributions are only / and /^cfacc- 
In the calculation of the diffuse components, the parameters 
and p^ are uniquely determined to fit the radio observa- 
tional data. The gas density nud, which affects the calculated 
flux of the gamma-ray emission, ranges within the values in 
Table O We note that the gamma-ray emission in the diffuse 
region is negligible within this range of nud, compared with 
the emission in the filamentary region. On the other hand, the 
gas properties of the filaments have free parameters of nuf and 
^atc- When we set nuf and n^tc within the values in Table[2l the 
parameters Bf, pi and / are uniquely determined to fit the 
observational data in the gamma-ray and radio bands. 

Figure [9] shows the spectra in the case of nm = 250 cm~^, 
^atc = 6 cm~^, and ^Hd = 0-4 cm~^. The cases of different 
values of the parameters are summarized in Table |5] Note 
that the spectra are almost the same for the values of the pa- 
rameters in the table. The gamma-ray flux is dominated by 
the TT^-decay emission from the dense filaments d ue to high 
densities of gas and CRs as in the case of S ection 14.11 This 
is consistent with the result in Section 13.11 which indicates 
the spatial correlation between the GeV gamma-ray and the 
filamentary Ha emissions. The obtained magnetic pressures 
are consistent with the values in Table [21 except for the case 
of nuf =100 cm~^ where Bf = 10 jiG is slightly smaller than 
that in the table. In all cases, the filaments are supported by 
thermal pressure. While the intensity of emission can be ex- 
plained well, the observed spectral index is harder than the 



calculated one at lower energy below the break in the ra- 
dio band. This discrepancy could be attributable to a cutoff 
around 100 MHz (4 x 10~^ eV) due to free-free absorption, 
which makes the spectral index harder than predicted by our 
model. 

In this paper, we assume the compression ratio of 4 based 
on the specific heat ratio of non-relativistic gas (=5/3). The 
compression ratio is expected to increase if we consider the 
effect of the relativistic gas whose specific heat ratio is 4/3. 
The precise treatment of this effect is beyond the scope of this 
paper. Here, let us consider the case of a compression ratio 
of 7 based on the specific heat ratio of the relativistic gas, i.e., 
CRs. In this case, the calculated spectral index of the particle 
energy distribution becomes harder by ~ 0.5. The reaccel- 
eration model can still explain the observed multiwavelength 
spectra when we set the values of / and to ~ 60% of those 
in the default case where the compression ratio is 4. The 
harder spectra explain the observed radio data better than the 
default case. 

Middle-aged SNRs interacting with molecular clouds con- 
stitute the dominant class o f SNRs detected by the Fermi 
LAT (s ee Uchiyama| 1201 1|) . It has been proposed by 
lUchiya ma et al. (201^ that the radiative filaments formed 
through interactions between molecular clouds and the SNR 
blast wave can account for the high gamma-ray luminosity of 
these SNRs (the Crushed Cloud model). In most cases, the 
reacceleration of Galactic CRs is sufficient to supply the re- 
quired CR density in the filaments. The scenario discussed in 
this section can be regarded as the atomic cloud version of the 
Crushed Cloud model, which is indeed insensitive to the phys- 
ical parameters of a pre- shock c loud. The GeV gam ma-ray 
emission from the Cygnus Loop ( iKatagiri et al.ll201 ih may be 
explained by the Crushed Cloud model as well. 




E [eV] 

Fig. 9. — SEDs of SNR S147 with model curves of emission from 
the reaccelerated CRs in the case of nm = 0.4 cm~^, nm = 250 cm~^, 
and natc = 6 cm~^. The distance to the SNR is assumed to be 
d = \3 kpc. The observed data in the gamma-ray and the radio 
bands (black points) are the same as Figure |7] The red and green 
lines represent diffuse and filamentary components, while the blue 
lines represent a sum of both fluxes. The radio emission is mod- 
eled by synchrotron radiation from the relativistic electrons (dashed 
line). The gamma-ray emission is explained with a combination of 
TT^-decay from the relativistic protons (dashed line), bremsstrahlung 
from primary (dot-dashed line) and secondary electrons (vertical- 
dashed line), and IC scattering from primary electrons (dotted line). 



4.3. Dependence on the Distance 

In Section |4~2l the distance to S147 is adopted as that to the 
plausibly associated pulsar. However, the absorption lines of 
the Ble star indicate a different distance. In this section, we 
consider a dependence on the dist ance for the reacceleration 
model as described in Section |4~2l 

We first evaluate the dependence on the distance of the 
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TABLES 

Parameters for the Model of Reaccelerated Cosmic-ray Emission (d = 13 kpc) 









Parameters 






Energetics 






[cm-3] 


[cm-3] 


B 
[/iG] 


Pc acc 

[GeV c-i] 


/ 

[%] 


[10^8 erg] 


We 

[10^6 erg] 


Up 
[eV cm-3] 


Ub 
[eV cm-3] 


filament 
diffuse 


100 
250 
500 
0.1-0.4 


2-6 
2-6 
2-6 


70 
120 
190 
13 


40-60 
30-40 
20-30 
8 


0.6-0.8 
0.2-0.3 
0.10-0.13 


0.36 
0.15 
0.080 
6.6 


1.8 
0.68 
0.30 

71 


80-270 
220-820 
480-1600 
2.5 


120 
360 
900 
4 



Note. — The total kinetic energies (£"1^111) of radiating protons (Wp) and electrons (We) are calculated for E'kin > 50 MeV and £'kin > 20 MeV, 
respectively. 

^ The densities of the gas are taken from Table|2] See Section |4?2| f or details. 



TABLE 6 

Parameters for the Model of Reaccelerated Cosmic-ray Emission (d = 0.88 kpc) 









Parameters 






Energetics 






nil 




B 


Pc acc 


/ 












[cm~^] 


[cm~^] 


[mG] 


[GeVc-1] 


[%] 


[10^8 erg] 


[10^6 erg] 


[eV cm ^] 


[eV cm~^] 


filament 


250 


10 


120 


50 


0.6 


0.070 


0.31 


100 


360 


dependence on d 












J2 


J2 






diffuse 


1.8 




17 


7 




1.8 


21 


2.2 


7 


dependence on d 












J3 


J3 







Note. — The total kinetic energies (£"1^111) of radiating protons (Wp) and electrons (Wg) are calculated for £'kin > 50 MeV and £'kin > 20 MeV, 
respectively. 



parameters in the diffuse region. The gas density in the 
diffuse region is estimated to be ^Hd = 4^icm oc using 
Equation |5l The total energy of the accelerated particles is 
W^^p (X- n'^eAP^y ^ The energy density = /V is 
thus independent of the distance. The emissivity of the syn- 
chrotron emission is proportional to assuming that 
the radiating electrons obey a power-law distribution with the 
index of m. Here we approximate the distributions of the 
accelerated particles by a power-law distribution with an in- 
dex of 1.8. The flux of the calculated synchrotron emission 
should explain the radio observational data for any distance; 

(the calculated flux) ex B^^'^^^^^W^/d'^ = constant. The mag- 
netic field is thus proportional to ~ d~^-^ . The photon 
energy of the synchrotron emission is approximately propor- 
tional to j^B. To explain the observational break in the radio 
band, oc^ 7 ex B'^^-^ cx J^-^^. Next, the parameters of 

the filaments are considered. The filamentary gas density ^Hf 
is derived from the optical observation, and independent of 
the distance. On the other hand, the gas density of the atomic 
cloud is ^atc d~^-^ using Equation [SI where we use an ap- 
proximate proportionality of F cx d^-^ around d=\3 kpc. The 
observed gamma-ray emission is reproduced by the tt^ -decay 
emission in the filaments; (the calculated flux) cx nm^l/d^ = 
constant, which leads the relation Wp oc d^. Since W^/d^ cx 
Wp/d'^ is independent of the distance, the magnetic field Bf 
should also be independent of the distance to reproduce the 
observed synchrotron emission in the radio band. Here we de- 
fine the compression ratio rcomp as ^Hf/^atc d^-^. The num- 
ber distribution of the accelerated particles has a dependence 
on distance of nlp(p) = s^^^n^cc(s~^^^p) 5'(2+'^)/3 j3.2^ 



where s = rcomp /^sh oc d^-^. The energy density is thus esti- 
mated to be Up cx nl p(p) cxd^-^. The total energy is described 
as cx /V/rcomp x ^l^p(p)^ which yields the filling factor 
of / oc d~^-^ . The cutoff momentum pi is determined by 
the observed gamma-ray data. Since pi is approximately 

proportional to the gamma-ray cutoff calculated by the tt^- 
decay emission, the cutoff has a dependence on distance of 

Table [6] shows one of the parameter sets of the reaccelera- 
tion model for the case of J = 0.88 kpc. The dependence on 
the distance d is also summarized in it. The observed inten- 
sity of the gamma-ray and radio emission is also reproduced 
by the reacceleration model for J = 0.88 kpc. We note that the 
contribution of gamma rays in the diffuse region is not negli- 
gible below a few GeV when ^Hd is greater than ~ 2 cm~^. In 
any case, the gamma-ray emission of the filaments is needed 
to explain the observed data above a few GeV. 

5. CONCLUSIONS 

We have presented results for the GeV gamma-ray observa- 
tions of the region around SNR SI 47 using about 31 months 
of data accumulated by the Fermi LAT. The gamma-ray emis- 
sion is spatially extended, being consistent with the size of 
S147 (~ 2000- There is no indication that the gamma-ray 
emission comes from the associated pulsar PSR 105 38+28 17. 
The best fit to the LAT data is obtained using the Ha im- 
age as a SNR spatial template rather than simple geometrical 
shapes. Comparisons between gamma-ray and Ha fluxes in- 
dicate a possible correlation between them, suggesting that 
the gamma rays come from the thin filaments observed in 
the Ha and radio bands. The observed energy spectrum be- 
tween 0.2-200 GeV has a indication of spectral steepening; a 
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smoothly broken power law provides a better fit than a sim- 
ple power law at 2cr significance, which means the possibility 
of a simple power law cannot be rejected. The gamma-ray 
luminosity amounts to 1.3 x 10^"^ (J/ 1.3 kpc)^ erg s"^ 

The LAT spectrum is best explained by tt^ -decay gamma 
rays from relativistic protons in the dense filaments. We find 
that the reacceleration of pre-existing CRs and subsequent 
adiabatic compression in the filaments is sufficient to provide 
the required energy density of high-energy electrons and pro- 
tons. There are two main distance estimates to SI 47: from 
the pulsar association (1.3 kpc) or from the absorption lines 
(< 0.88 kpc). We consider the cases of 1.3 kpc and 0.88 kpc, 
and the gamma-ray emission can be explained in the same 
framework for either distance. SNR S147 offers a firm exam- 
ple of the realization of the Crushed Cloud model, and sup- 
ports the importance of dense filaments in SNRs as gamma- 
ray production sites. 
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